Int J Life Cycle Assess (2014) 19:1469-1481 
DOI 10.1007/s 11367-014-0755-9 


CARBON FOOTPRINTING 


An appraisal of carbon footprint of milk from commercial 
grass-based dairy farms in Ireland according to a certified life 
cycle assessment methodology 

Donal O’Brien • Padraig Brennan • James Humphreys • 

Eimear Ruane • Laurence Shalloo 


Received: 3 February 2014 / Accepted: 11 May 2014 /Published online: 14 June 2014 
© The Author(s) 2014. This article is published with open access at Springerlink.com 


Abstract 

Purpose Life cycle assessment (LCA) studies of carbon foot¬ 
print (CF) of milk from grass-based farms are usually limited 
to small numbers of farms (<30) and rarely certified to inter¬ 
national standards, e.g. British Standards Institute publicly 
available specification 2050 (PAS 2050). The goals of this 
study were to quantify CF of milk from a large sample of 
grass-based farms using an accredited PAS 2050 method and 
to assess the relationships between farm characteristics and 
CF of milk. 

Materials and methods Data was collected annually using on- 
farm surveys, milk processor records and national livestock 
databases for 171 grass-based Irish dairy farms with informa¬ 
tion successfully obtained electronically from 124 farms and 
fed into a cradle to farm-gate LCA model. Greenhouse gas 
(GHG) emissions were estimated with the LCA model in C0 2 
equivalents (C0 2 -eq) and allocated economically between 
dairy farm products, except exported crops. Carbon footprint 
of milk was estimated by expressing GHG emissions attribut¬ 
ed to milk per kilogram of fat and protein-corrected milk 
(FPCM). The Carbon Trust tested the LCA model for non¬ 
conformities with PAS 2050. PAS 2050 certification was 
achieved when non-conformities were fixed or where the 
effect of all unresolved non-conformities on CF of milk was 
<±5 %. 

Results and discussion The combined effect of LCA model 
non-conformities with PAS 2050 on CF of milk was <1 %. 
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Consequently, PAS 2050 accreditation was granted. The mean 
certified CF of milk from grass-based farms was 1.11 kg of 
C0 2 -eq/kg of FPCM, but varied from 0.87 to 1.72 kg of C0 2 - 
eq/kg of FPCM. Although some farm attributes had stronger 
relationships with CF of milk than the others, no attribute 
accounted for the majority of variation between farms. How¬ 
ever, CF of milk could be reasonably predicted using N 
efficiency, the length of the grazing season, milk yield/cow 
and annual replacement rate (R =0.75). Management changes 
can be applied simultaneously to improve each of these traits. 
Thus, grass-based farmers can potentially significantly reduce 
CF of milk. 

Conclusions The certification of an LCA model to PAS 2050 
standards for grass-based dairy farms provides a verifiable 
approach to quantify CF of milk at a farm or national level. 
The application of the certified model highlighted a wide 
range between the CF of milk of commercial farms. However, 
differences between farms’ CF of milk were explained by 
variation in various aspects of farm performance. This implies 
that improving farm efficiency can mitigate CF of milk. 

Keywords Carbon footprint • Certification • Dairy • Grass • 
Greenhouse gas - Milk 

1 Introduction 

Grass-based dairy production is a key agricultural industry in 
some developed nations, particularly Ireland and New 
Zealand. The dairy sector internationally, however, is also an 
important source of greenhouse gas (GHG) emissions, respon¬ 
sible for approximately 3 % of global emissions (Opio et al. 
2013), 10 % of Ireland’s emissions and up to 20 % of New 
Zealand emissions (Beukes et al. 2010; Deighton et al. 2010). 
Climate change caused by GHG emissions has become an 
important political issue. This has led to growing awareness 
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amongst consumers of the potential adverse effects of climate 
change, which is expected to increase the demand for food 
products that generate low GHG emissions (Roy et al. 2009). 
In light of this anticipated demand, major retailers (e.g. Tesco 
and Wal-Mart) have already put in place sustainability pro¬ 
grams to monitor emissions associated with the production of 
their food products. Therefore, quantifying GHG emissions 
from all milk-producing nations is becoming a pre-requisite, 
but especially for countries like Ireland, which export the 
majority (85 %) of their dairy products (CSO 2013). 

Simple simulation models of complex biological and tech¬ 
nical processes are used to quantify GHG emissions from milk 
given that direct measurement of agricultural GHG sources 
(e.g. soils) is difficult and cost prohibitive. The preferred 
approach to simulate GHG emissions from milk production 
is the life cycle assessment (LCA). Ideally, LCA should be 
applied to quantify emissions from all life cycle phases of 
milk, such as farming, processing, consumption and waste 
treatment. The methodology, though, is generally applied to 
the cradle to farm-gate stage of dairy production and, thus, 
only simulates emissions associated with milk prior to the sale 
of the product from the farm (Yan et al. 2011). Firstly, this is 
because of the variation between farming systems, and sec¬ 
ondly, most GHG emissions from milk are emitted during the 
farm stage. For instance, LCA studies of milk to the retail 
stage report the cradle to farm-gate stage causes 78-95 % of 
GHG emissions from milk (Defra 2007; Gerber et al. 2010). 

Usually studies apply a cradle to farm-gate LCA approach 
to compare the GHG emission intensity or carbon footprint 
(CF) of milk (kg of GHG/unit of milk) from contrasting dairy 
farms (O’Brien et al. 2010; Flysjo et al. 2011) or use the 
approach to assess GHG mitigation strategies (Rotz et al. 
2010). The majority of previous research though has been 
carried out on a limited number (<20-3 0) of farms (generally 
not randomly chosen) or based on a single-point modelled 
national average farm, because significant resources are re¬ 
quired to conduct an LCA study of dairy farms (Thomassen 
et al. 2009). However, researchers are increasingly overcom¬ 
ing this resource challenge, in part through improvements in 
mobile computer technology (e.g. Thoma et al. 2013). Thus, 
this implies that there is an increasing scope to complete LCA 
studies of CF of milk with a large sample of farms (>100). 

The major advantage of performing LCA on a large num¬ 
ber of dairy farms is that it provides an insight into the 
variation between CF of milk from commercial farms, which 
is generally not possible through national average estimates or 
small-scale studies. In addition, differences between the CF of 
milk from dairy farms may be related to variation in farm 
performance (e.g. milk yield/ha), which can facilitate the 
development of management practices or mitigation strategies 
to reduce CF of milk (Thomassen et al. 2009). However, apart 
from DairyCo (2012), previous large-scale European cradle to 
farm-gate LCA studies of CF of milk have only considered 


farms where cows do not graze or graze for a short period, 
usually no longer than 5 months. Thus, strategies suggested to 
reduce CF of milk from such studies will not be applicable to 
grass-based farms where cows graze for an extended period 
(e.g. 9-10 months; O’Brien et al. 2012). 

Although several studies have used LCA to quantify CF of 
milk from grass-based dairy farms, the application of the 
methodology varies between studies (Yan et al. 2011). To try 
to overcome this methodology challenge, LCA was applied 
according to the British Standards Institute (BSI 2011) pub¬ 
licly available specification 2050:2011 (PAS 2050) LCA stan¬ 
dard for GHG emissions. In addition, to comply with PAS 
2050 certification requirements, an accredited third party 
(Carbon Trust) assessed all LCA procedures. The goal of this 
study was to quantify CF of milk from a significant group 
(>100) of grass-based dairy farms using an LCA method 
independently accredited to comply with PAS 2050 standards. 
Additionally, we aimed to determine the factors that caused 
variation amongst farms’ carbon footprints of milk. The study 
was limited to Ireland, but is also relevant to developed 
nations where extended grazing of dairy cattle is practiced, 
given that the application of a certified LCA method to assess 
emissions from a large sample of grass-based dairy farms is 
rare. 

2 Materials and methods 

2.1 Data collection 

A random sample of 171 commercial Irish dairy farms was 
audited between November 2011 and December 2012. The 
farms sampled were located in the northeast, east, southeast, 
south and southwest of the country. The farms assessed were 
therefore a subset of the country’s dairy farms and, thus, not 
representative of the population’s CF of milk. The Carbon 
Trust provided guidance on randomly selecting dairy farms 
for the region. All data was collected electronically using 
annual on-farm surveys, electronic feeds of dairy processor 
milk data and livestock data (DAFM 2011; ICBF 2012). 

The Department of Agriculture, Food and the Marine 
(DAFM) animal identification and movement (AIM) system 
was used to source bovine data from dairy farms. The AIM 
system records all births, movements and disposal of farm 
animals on a daily basis. It is highly reliable and used to meet 
EU requirements on traceability for bo vines (DAFM 2011). 
Trained auditors carried out farm surveys that corresponded to 
a 1-year period (fiscal year). Three to 4 hrs of data collection 
was sufficient to survey a farm, provided source documents, 
for instance accounting books, were available. Farm auditors 
collected information on key parameters, for instance, the area 
of the dairy enterprise, grazing management, concentrate feed, 
manure management and use of fertiliser, fuel and contractors. 
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In addition, auditors collected farm livestock inventory data, 
which was cross-checked with national livestock databases 
(DAFM 2011; ICBF 2012). 

Electronic data collected was fed into a farm database to 
operate an LCA model to calculate the CF of milk from dairy 
farms. The process described to collect data was externally 
verified by the Carbon Trust. The Carbon Trust audited farms 
by cross-checking electronic data collected with farm in¬ 
voices, milk supply records, auditor measurements and farm 
livestock records. During this validation procedure, 45 farms 
were excluded from the analysis because of inadequate or 
unreliable data (e.g. inconsistency between farmer records of 
livestock numbers and AIM records). Two more farms were 
omitted from the evaluation because they had or began to 
cease producing milk during the period of the study. In total, 
124 dairy farms were analysed. 

The majority (70 %) of farms analysed were grass-based 
spring calving systems where the aim was to minimize costs 
through maximizing the proportion of grazed grass in the diet 
(Kennedy et al. 2005). To achieve this goal, farmers synchro¬ 
nized calving with the onset of grass growth in early or mid¬ 
spring and cows remained at pasture from calving until late 
autumn or early winter. Pasture was usually offered to cows 
through a rotational grazing system, where cows were offered 
sections of pasture for 1-2 days or until a specific grazing 
height was reached (e.g. 4-5 cm) and then moved to a new 
section. When grass growth exceeded herd feed demand, 
surplus grass was harvested as grass silage, hay or both and 
fed to cows indoors from early winter to early spring. Cows 
were supplemented with purchased concentrate feeds, when 
grass growth was insufficient to meet herd feed requirements. 

Spring calving grass-based dairy systems are the dominant 
method of producing milk in Ireland (Evans et al. 2004). Thus, 
the supply pattern of milk is highly seasonal, greatest in May 
and lowest in January. There is a requirement for some milk to 
be produced out of season for the fluid milk market and 
specialty-type products. In order to incentivize farmers to 
supply milk during the winter period, dairy processors offer 
farmers a winter contract with a milk price bonus to cover the 
extra costs of production (Fitzgerald et al. 2004). To meet this 
requirement, a minority of farms analysed calved cows in 
autumn or throughout the year. Thus, these farms offered 
milking cows concentrate feed over the winter and, in addition 
to grass silage, fed maize silage, whole crop silage or both. 
Farms that milked cows throughout the winter, however, also 
aim to graze cows for an extended period (6-8 months) to 
minimize costs, particularly in late autumn and early spring 
(Fitzgerald et al. 2004). 

2.2 Greenhouse gas simulation 

A cradle to farm-gate attributional LCA model developed by 
O’Brien et al. (2010) was used to simulate annual GHG 


emissions from dairy farms. Thus, all sources of GHG emis¬ 
sions associated with dairy production until milk leaves the 
farm were simulated, including off-farm GHG sources such as 
fertiliser production (Fig. 1). The model calculated GHG 
emissions by combining information from the farm database, 
mentioned previously, with literature emission algorithms. 
However, the LCA model procedures and literature emission 
factors were adapted based on farm data availability and 
recommendations of the Carbon Trust. In addition, fluorinated 
gases (F-gases) from refrigerant loss were added to the model 
as a source of GHG emissions and estimated based on farm 
service records of cooling equipment. 

The main on-farm sources of GHG emissions quantified by 
the LCA model are summarized in Table 1 . Enteric methane 
(CH 4 ) emissions were estimated by firstly computing animal 
feed intakes. To ensure animal feed intakes were realistic, they 
were calculated to fulfil net energy requirements for milk 
production, maintenance, pregnancy and body weight change 
(Jarrige 1989). Where possible, information directly acquired 
on animal and feed variables were used to validate feed 
intakes. However, this was not possible for all variables, for 
instance energy values of forages. Thus, data from literature 
sources were also used (O’Mara 1996, 2006). Enteric CH 4 
emissions were estimated as a function of intake, using a fixed 
factor of 6.5 % of gross energy intake (GEI) when cattle 
grazed grass (Duffy et al. 2011b). However, when the diet 
comprised of only silage and concentrate, a regression 
equation from Yan et al. (2000) was used. 

Methane emissions from manure were estimated according 
to the Intergovernmental Panel on Climate Change IPCC 
(2006) guidelines as a proportion of the maximum CH 4 po¬ 
tential ( B 0 ) of manure volatile solids (VS). Based on O’Mara 
(2006), manure VS excretion was calculated by multiplying 
animal organic matter (OM) intake by the indigestible OM 
component of the diet. The quantity of manure VS requiring 
storage was computed based on the number of days animals 
spent housed. The B 0 of manure VS for dairy cattle was based 
on Duffy et al. (2011b). The proportion of the B 0 that was 
emitted from manure VS was computed using specific CH 4 
conversion factors for different manure storage systems, but 
for manure deposited on pasture, a default factor (1 %) was 
used (IPCC 2006). The CH 4 conversion factors for manure 
storage systems were selected assuming an annual average 
ambient temperature of 10 °C for Ireland (Met 2013). 

Emissions of nitrous oxide (N 2 0) from manure were de¬ 
rived after calculating N excretion, which was estimated as the 
difference between total N intake and N output in meat and 
milk. Direct N 2 0 emissions from manure storage were esti¬ 
mated using manure storage-specific emission factors. After 
subtraction of N losses during housing and storage, N 2 0 
emissions from manure spreading were estimated as 1 % of 
N applied (IPCC 2006). In addition, this factor was used to 
estimate N 2 0 emissions from synthetic fertiliser spreading, 
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Fig-1 An illustration of the 
major sources of on- and off-farm 
greenhouse gas emissions, carbon 
dioxide (C0 2 ), methane (CH 4 ), 
nitrous oxide (N 2 0) and 
fluorinated gases (F-gases), 
quantified using a cradle to farm- 
gate life cycle assessment model 
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crop residues and from soil mineralisation following land use 
change. Nitrogen inputs from crop residues and soil 
mineralisation were computed based on the IPCC (2006) 
guidelines. A greater N 2 0 emission factor was estimated for 
manure excreted by grazing cattle (2 % of N excreted), com¬ 
pared to manure or fertiliser spreading, given that urine de¬ 
posited by grazing cattle can cause large N losses (Van 
Groenigen et al. 2005). 

Indirect N 2 0 emissions from re-deposition of volatilized 
NH 3 were computed as 1 % of NH 3 emitted from fertiliser and 
manure (IPCC 2006). Volatilization of NH 3 from synthetic 
fertiliser application was estimated using N loss factors from 
Hyde et al. (2003). A mass flow approach and emission 
factors from Hyde et al. (2003) were used to estimate NH 3 
emitted during cattle housing, manure storage and spreading. 
Ammonia loss from manure excreted by grazing cattle was 
estimated as 20 % of N excreted on pasture (IPCC 2006). 
Indirect N 2 0 emitted from leaching of N was estimated as 
0.75 % of N leached (IPCC 2006). A default leaching factor 
from Duffy et al. (2011b) was used to estimate nitrate 
leaching. Nitrogen available for leaching was quantified by 
subtracting N removed in products and N lost directly to the 
atmosphere from total N inputs (manure, soil mineralisation, 
fertiliser and crop residues). 

On-farm emissions of carbon dioxide (C0 2 ) from fossil 
fuels, lime and urea were estimated using the IPCC (2006) 
guidelines. Short-term biogenic sources and sinks of C0 2 
such as animals, crops and manure were considered to be 


neutral with respect to GHG emissions given that the IPCC 
(2006) assume all C absorbed by animals, crops and manure to 
be quickly released back to the atmosphere through respira¬ 
tion, burning and decomposition. Agricultural soils also have 
the potential to emit or sequester C0 2 (Rotz et al. 2010). 
However, to comply with PAS 2050, C sequestration was 
not included for permanent pasture. This was because the 
standard follows the IPCC (2006) recommendation that soil’s 
ability to store or lose C reaches equilibrium after a fixed 
period (20 years). Thus, land use change emissions were also 
restricted to this period and estimated as 6.7-7.0 t C0 2 /ha per 
annum when permanent on-farm grassland was converted to 
cropland (Carbon Trust 2013). 

Data acquired on external farm inputs (e.g. diesel and 
pesticides) were used primarily with emission factors (Table 2) 
from the Carbon Trust (2013) to estimate off-farm GHG 
emissions. For electricity generation and some other sources 
though, it was more appropriate to use emission algorithms 
from national literature sources. When emissions from an 
external farm input could not be estimated via the Carbon 
Trust or national literature, Ecoinvent (2010) was used. The 
production of specific imported feeds, for instance Malaysian 
palm kernel, was estimated to cause land use change emis¬ 
sions by computing the average land use change emissions for 
that crop in that country (Carbon Trust 2013). Directly attrib¬ 
uting land use change emissions to a crop conforms to the 
method used by the Food and Agriculture Organization (FAO) 
to estimate GHG emissions from milk (Opio et al. 2013). 


<£) Springer 












Int J Life Cycle Assess (2014) 19:1469-1481 1473 


Table 1 Key emission factors applied in a 
ammonia emissions and nitrate loss 

cradle to farm-gate life cycle assessment model to quantify on-farm greenhouse gas (GHG) emissions, 

Emission and source 

Emission factor 

Unit 

References 

Carbon dioxide (C0 2 ) 

Lime 

0.12 x lime application 

kg/kg lime 

IPCC (2006) 

Urea 

0.20 xurea application 

kg/kg urea 

IPCC (2006) 

Diesel 

2.63 xdiesel use 

kg/l 

IPCC (2006) 

Gasoline 

2.30xgasoline use 

kg/1 

IPCC (2006) 

Kerosene 

2.52 xkerosene use 

kg/l 

IPCC (2006) 

Liquefied petroleum gas (LPG) 

1.49 x LPG use 

kg/l 

IPCC (2006) 

Methane (CH 4 ) 

Enteric fermentation 

Daily cow and heifer (housing) 

DEIx (0.096+0.035 xy DMI /T DMI )-(2.298 xFL-1) 

MJ/day 

Yan et al. (2000) 

Daily cow and heifer (grazing) 

0.065 x GEI 

MJ/day 

Duffy et al. (2011b) 

Manure storage and excretion on pasture 

Manure VS excreted x 0.24x 0.67 x MS a x MCF b 

kg/year 

IPCC (2006); Met (2013) 

Nitrous oxide (N 2 0) 

Slurry storage 

0.005 x slurry N stored 

kg/kg N 

IPCC (2006) 

Solid manure storage 

0.005 x solid manure N stored 

kg/kg N 

IPCC (2006) 

Manure excreted on pasture 

0.02 xN excreted on pasture 

kg/kg N 

IPCC (2006) 

Synthetic N fertiliser 

0.01 X N fertiliser 

kg/kg N 

IPCC (2006) 

Manure application 

0.01 x (manure N applied-N storage loss) 

kg/kg N 

IPCC (2006) 

Crop residues 

0.01 xN crop residues 

kg/kg N 

IPCC (2006) 

Nitrate leaching 

0.0075 xN leached 

kg/kg N0 3 -N 

IPCC (2006) 

Ammonia (NH 3 ) re-deposition 

0.01 xsum of NH 3 loss 

kg/kg NH 3 -N 

IPCC (2006) 

Ammonia (NH 3 -N) 

Housing 

11—38 c 

g/luuk per day 

Duffy et al. (2011a) 

Slurry storage 

(2—4) d x area slurry store 

g/m 2 per day 

Duffy et al. (2011a) 

Solid manure storage 

94 x area solid manure store 

g/m 2 per year 

Duffy et al. (2011a) 

Slurry application 

(0.15-0.59) e xTAN in slurry spread 

kg/kg TAN 

Duffy et al. (2011a) 

Solid manure application 

0.81 xTAN in solid manure spread 

kg/kg TAN 

Duffy et al. (2011a) 

Grazing cattle 

0.2 x N excreted on pasture 

g/luuk per day 

Duffy et al. (2011a) 

Synthetic N fertiliser 

(0.08-0.23) f xN fertiliser applied 

kg/kg N 

Duffy et al. (2011a) 

Nitrate (N0 3 -N) 

N leaching 

0.1 x(N applied-NH 3 loss-N 2 0 loss) 

kg/kg N 

Duffy et al. (2011b) 


DEI digestible energy intake, S DMI silage diy matter intake, T DMI total diy matter intake, FL feeding levels above maintenance energy requirement, GEI 
gross energy intake, VS volatile solids, luuk UK livestock unit (equivalent to 500 kg body weight), TAN total ammoniacal nitrogen 

a MS = percentage of manure volatile solids managed in a specific storage system or percentage of manure excreted on pasture 

b MCF = methane conversion factor for manure volatile solids managed in a particular storage system or excreted on pasture. The MCF values were 0.01 
for manure excreted on pasture, 0.02 for solid manure system (dry matter (DM) >20 %), 0.66 for lagoon system and 0.17 for slurry system without a 
surface crust or 0.11 with a surface crust 

c Dependent on the age of the animal 

d Dependent on manure storage facility 

e Dependent on DM of slurry and season of application 

F 

Dependent on fertiliser compound 


2.3 Quantification and accreditation of CF of milk from dairy 
farms 

On- and off-farm GHG emissions were converted to C0 2 
equivalent (C0 2 -eq) emissions using the IPCC (2007) guide¬ 
lines’ global warming potential (GWP) factors, which have 


been revised (IPCC 2013) and summed to compute dairy 
farms’ annual C0 2 -eq emissions. The GWP factors for key 
GHG emissions were 1 for C0 2 ,25 for CH 4 and 298 for N 2 0, 
assuming a 100-year time horizon. The CF of milk from dairy 
farming was estimated by firstly allocating GHG emissions 
between farm outputs, milk, crops, manure and meat from 
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Table 2 Key emission factors used in a cradle to farm-gate life cycle assessment model for quantification of off-farm greenhouse gas emissions in 
kilograms of C0 2 equivalent 

Item 

Emission factor 

Reference 

Electricity, kWh 

0.60 

Howley et al. (2011) 

Diesel, 1 

0.41 

Ecoinvent (2010) 

Gasoline, 1 

0.57 

Ecoinvent (2010) 

Kerosene, 1 

0.39 

Ecoinvent (2010) 

Liquefied petroleum gas, 1 

0.30 

Ecoinvent (2010) 

Ammonium nitrate, kg N 

7.11 

Carbon Tmst (2013) 

Urea, kg N 

3.07 

Ecoinvent (2010) 

Lime, kg 

0.15 

Carbon Tmst (2013) 

P fertiliser, kg P 2 0 5 

1.86 

Carbon Tmst (2013) 

K fertiliser, kg K 2 0 

1.77 

Carbon Tmst (2013) 

Refrigerant, kg 

11.00-393.00 

Little (2002) 

Detergent, kg active ingredient 

0.11-1.03 

Ecoinvent (2010) 

Pesticide, kg active ingredient 

7.37 

Carbon Tmst (2013) 

Barley, kg dry matter (DM) 

0.35 

Carbon Tmst (2013) 

Com grain, kg DM 

0.45 

Carbon Tmst (2013) 

Citrus pulp, kg DM 

0.06 

Carbon Tmst (2013) 

Com gluten, kg DM 

0.34 

Carbon Tmst (2013) 

Molasses, kg DM 

0.15 

Carbon Tmst (2013) 

Rapeseed meal, kg DM 

0.40 

Carbon Tmst (2013) 

South America soybean meal, kg DM 

11.65 

Carbon Tmst (2013) 

USA soybean meal, kg DM 

0.32 

Carbon Tmst (2013) 

South America soybean hulls, kg DM 

0.28 

Ecoinvent (2010) 
Carbon Tmst (2013) 

USA soybean hulls, kg DM 

0.01 

Ecoinvent (2010) 
Carbon Tmst (2013) 

Compound concentrate, 16 % cmde protein a (CP), kg DM 

0.34 

Ecoinvent (2010) 
Carbon Tmst (2013) 

Compound concentrate, 20 % CP b , kg DM 

1.98 

Ecoinvent (2010) 
Carbon Tmst (2013) 

Compound concentrate, 32 % CP c , kg DM 

3.37 

Ecoinvent (2010) 
Carbon Tmst (2013) 


a Concentrate formulation on a DM basis: USA soy hulls 18 %, USA dried distillers grains 17 %, USA citrus pulp 17 %, USA com gluten feed 8 %, 
French rapeseed meal 8 %, German com grain 6 %, Cuban molasses 5 %, Irish barley 5 %, Malaysian palm kernel meal 4 %, Irish wheat feed 4 %, 
vegetable oil 3 %, Irish lime 3 %, minerals and vitamins 2 % 

b Concentrate formulation on a DM basis: Irish wheat feed 17 %, USA soy hulls 16 %, French rapeseed meal 15 %, Brazilian soybean meal 14 %, 
German com grain 7 %, Cuban molasses 6 %, Malaysian palm kernel meal 6 %, USA dried distillers grains 4 %, USA citrus pulp 4 %, French sunflower 
meal 4 %, Irish lime 3 %, vegetable oil 2 %, minerals and vitamins 2 % 

c Concentrate formulation on a DM basis: Brazilian soybean meal 26 %, French rapeseed meal 24 %, USA dried distillers grains 12 %, Malaysian palm 
kernel meal 8 %, USA soy hulls 7 %, Cuban molasses 6 %, USA com gluten feed 5 %, French sunflower meal 4 %, Irish lime 3 %, vegetable oil 3 %, 
minerals and vitamins 2 % 


culled cows and surplus calves. Where possible, PAS 2050 
recommendations to avoid allocation of GHG emissions be¬ 
tween products were applied. However, this was only 
achieved for exported crops by constraining the LCA model 
to quantify emissions from crops grown for dairy cattle. 

When allocation was required, GHG emissions were allo¬ 
cated between dairy farm products based on their economic 
value. The economic method of allocation was used instead of 
alternative methods, e.g. the physical allocation approach 


recommended by the IDF (2010), because economic alloca¬ 
tion is preferred by PAS 2050 when allocation cannot be 
avoided. The economic value of milk and meat was estimated 
using 5-year average prices from 2007 to 2011 (CSO 2013). 
Cattle manure exported from dairy farms had no economic 
value. Thus, GHG emissions were allocated to cattle manure 
based on the site of storage and application. Emissions from 
the transport of exported manure were attributed to the im¬ 
porter. Greenhouse gas emissions attributed to milk were 
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expressed per kilogram of fat and protein-corrected milk 
(FPCM) to quantify the CF of milk from dairy farms. An 
algorithm from the International Dairy Federation (IDF 2010) 
was used to estimate FPCM as 4 % fat and 3.3 % tme protein. 

The LCA model estimate of CF of milk was tested by the 
Carbon Trust to certify compliance with PAS 2050. Research 
data from Teagasc (2011) on poor, average and high 
performing grass-based Irish dairy farms were initially used 
to evaluate the model. Subsequently, research data from 
Olmos et al. (2009) was used to test the model’s ability to 
estimate CF of milk from a high input farm where cows did 
not graze. The LCA model was also tested with data from 
commercial farms. This was achieved by randomly selecting 
12 grass-based farms from the study sample. 

To achieve accreditation, non-conformities identified by 
the Carbon Trust between the LCA model and PAS 2050 were 
addressed. Non-conformities that had a non-material impact 
on the CF of milk were either justified through research data 
or changed based on advice of the Carbon Trust. PAS 2050 
certification was granted when all non-conformities were 
addressed or where the effect of all unresolved non-material 
non-conformities on the CF of milk was<±5 %. To determine 
if the LCA model CF result was within the<±5 % threshold, 
all changes recommended by the Carbon Trust were applied to 
any model non-conformity with PAS 2050. 

2.4 Statistical analyses 

The Statistical Analysis Systems (SAS) Institute software 
package (SAS 2008) was used to evaluate relationships be¬ 
tween farm characteristics (e.g. herd size) and CF of milk. The 
strength of the relationship between CF of milk and individual 
farm characteristics was measured using the Pearson correla¬ 
tion where the data was normally distributed and the Spear¬ 
man rho correlation for non-normal data. Correlations were 
also performed amongst farm characteristics, where a farm 
variable was correlated to CF of milk. Least squares regres¬ 
sion analysis was used to evaluate the associations between 
individual farm variables and CF of milk. The equality of the 
variances of residuals of regression models was checked for 
normality visually and the Shapiro-Wilk test was used 
(^<0.05). The significance of regression coefficients was 
assessed using the t statistic. 

Non-linear terms were added to regression models based 
on a visual assessment of the distribution of the residuals. The 
stepwise multiple regression procedure of SAS (2008) was 
used to determine whether the addition of non-linear terms 
made a significant contribution (,P<0.05) to a regression 
equation. The procedure was also used to develop predictive 
models by assessing the relationship between the CF of milk 
and all farm characteristics listed in Tables 3 and 4. Variance 
inflation factors and condition indices were applied in SAS 
(2008) to check for multicollinearity within multiple 


regression models. Regression coefficients with variance in¬ 
flation factors >10 were omitted. 


3 Results 

3.1 General farm characteristics 

The weighted mean, standard deviation (SD), coefficient of 
variation (CV) and range for various farm characteristics of 
124 dairy farms are shown in Table 3. On average, dairy farms 
were 51 ha in size, stocked at 2.1 livestock units (LU)/ha and 
produced 489 t of FPCM. The mean herd size was 91 cows 
and the replacement rate averaged 19 % (SD=8 %) across 
farms. Cows spent the majority of the year at pasture (mean of 
245 grazing days, SD=23). Thus, their diet was mainly com¬ 
posed of grazed grass (mean 63 %; SD=7 %). Concentrate 
input averaged 643 kg DM/cow (SD=310) across farms, and 
the mean quantity of forage purchased was 118 kg DM/cow 
(SD = 202). On-farm N fertiliser application averaged 
172 kg N/ha (SD=54) and the mean farm-gate N surplus (N 
imports-N exports) was 150 kg N/ha (SD=53). On-farm 
electricity consumption averaged 47 kWh of electricity/t of 
FPCM (SD=18) and fuel use was 97 1 of fuePha (SD=15). 

3.2 Certified CF of milk from Irish dairy farms 

All non-conformities identified between the LCA model and 
PAS 2050 during evaluation of the CF of milk from dairy 
farms were addressed. Generally, non-conformities were re¬ 
solved by applying changes recommended by the Carbon 
Trust, but for some non-conformities considered non¬ 
material (e.g. simulation of manure excretion by cattle), no 
revisions were undertaken. Overall, the cumulative effect of 
non-material non-conformities was <1 %. Thus, CF of milk 
estimated by the LCA model was accredited to comply with 
PAS 2050. 

On average, 90 % of dairy farms’ GHG emissions were 
allocated to milk production. The mean certified CF of milk to 
the farm-gate, weighted by farm milk production, was 1.11 kg 
of C0 2 -eq/kg of FPCM. The SD was within ±0.13 kg of C0 2 - 
eq of the weighted mean. Figure 2 shows across the 124 dairy 
farms that there was significant variability between farms’ CF 
of milk, ranging from 0.87 to 1.72 kg of C0 2 -eq/kg of FPCM. 
On average, CF of milk mainly consisted of CH 4 (47 %) and 
N 2 0 (34 %) emissions, followed by emissions of C0 2 (19 %) 
and F-gases (<0.5 %). Approximately, 80 % of the CF of milk 
to the farm-gate was caused by GHG emissions generated 
directly on-farm. 

The largest individual on-farm contributor to CF of milk 
was CH 4 from enteric fermentation (44 %), which averaged 
113 kg CH 4 /cow. Nitrous oxide from manure excreted by 
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Table 3 Weighted means, standard deviations (SD), coefficients of variations (CV), minimum, maximum and lower and upper 10 percentiles of key 
farm characteristics of 124 Irish dairy farms 


Farm characteristic 

Mean 

SD 

CV (%) 

Min 

Lower 10 % 

Upper 10 % 

Max 

Farm size, ha 

51 

27 

52 

18 

30 

75 

245 

Stocking rate, LU/ha 

2.10 

0.46 

22 

1.30 

1.50 

2.60 

2.93 

Cows, n 

91 

55 

60 

29 

44 

149 

468 

Annual FPCM yield/cow 

5,380 

834 

16 

3,163 

4,433 

6,483 

8,071 

Annual FPCM yield/ha 

9,494 

2,534 

27 

4,529 

6,435 

12,832 

15,085 

Fat, % 

3.89 

0.17 

4 

3.54 

3.71 

4.13 

4.31 

Protein, % 

3.38 

0.10 

3 

3.10 

3.26 

3.51 

3.73 

Lactose, % 

4.70 

0.05 

1 

4.48 

4.64 

4.76 

4.84 

Replacement rate, % 

19 

8 

43 

7 

9 

28 

36 

Grazing season, days 

245 

23 

9 

184 

214 

274 

280 

Concentrate, kg DM/cow 

643 

310 

48 

224 

310 

1,062 

1,552 

Purchased forage, kg DM/cow 

118 

202 

171 

0 

0 

337 

1,197 

Total feed, kg DM/cow 

5,630 

549 

10 

4,270 

4,961 

6,350 

6,444 

On-farm N fertiliser, kg/ha 

172 

54 

32 

95 

112 

255 

307 

Farm-gate N surplus a , kg/ha 

150 

53 

36 

70 

87 

231 

240 

Farm N efficiency 13 , % 

27 

7 

25 

14 

20 

37 

39 

Electricity, kWh/cow 

253 

98 

38 

57 

152 

388 

409 

On-farm fuel, 1/ha 

97 

15 

15 

41 

78 

116 

121 


LU livestock unit (equivalent to the average annual N excretion of an Irish daily cow); FPCM fat and protein-corrected milk, where milk was 
standardised to 4 % fat and 3.3 % true protein per kilogram; DM dry matter 

a N imports-N exports passing in or out through the farm-gate 

b Farm N exports/farm N imports 

grazing cattle was the next largest on-farm contributor (17 %) 
to CF of milk, followed by N 2 0 and C0 2 emissions from 
fertiliser application (8 %), CH 4 and N 2 0 emissions from 
manure storage and spreading (6 %) and C0 2 emissions from 
fuel consumption and lime (5 %). The remainder of the CF of 
milk comprised of off-farm GHG emissions from fertiliser 
manufacture (10 %), concentrate production (7 %) and energy 
generation (3 %). The majority of GHG emissions from 
concentrate production were caused by C0 2 emissions from 
land use change (57 %). 

3.3 Associations between farm characteristics and CF of milk 

Significant (^<0.05) associations between farm characteris¬ 
tics and CF of milk were linear. Farm N efficiency had the 
strongest negative correlation with CF of milk (r=-0.66; 

Table 4), which regression analysis found decreased CF of 
milk by 13 g of C0 2 -eq (standard error (se)= 1) for each 1 % 
increase in N efficiency. Apart from total farm milk produc¬ 
tion, measures of milk yield were moderately negatively cor¬ 
related to CF of milk, but the correlation was stronger when 
evaluated per hectare (r=-0.55) rather than per cow (r= 

-0.48). Milk yield/cow had a similar negative correlation with 
CF of milk as homegrown forage utilisation per hectare, but a 
slightly stronger correlation than the length of the grazing 


season (r=-0.45). Moderate to weak negative correlations 
(r=-0.24 to -0.32) occurred between CF of milk and milk 
fat content, milk protein content, total feed intake/cow and 
stocking rate. 

Farm N surplus per hectare and purchased forage/cow were 
the only farm attributes positively correlated to CF of milk, but 
the correlations were weak (r=0.20 to 0.28). Annual replace¬ 
ment rate, farm size and herd size were not associated with CF 
of milk. However, stepwise multiple regression (Table 5) 
showed that CF of milk could be reasonably well explained 

'-y 

(7C=0.75) by farm N efficiency, the length of the grazing 
season, milk yield/cow and annual replacement rate. Of these 
farm attributes, the procedure found that the majority of var¬ 
iation in CF of milk amongst farms was explained by farm N 
efficiency. 


4 Discussion 

The application of LCA according to PAS 2050 highlighted a 
significant variation in the CF of milk amongst grass-based 
dairy farms. The differences between CF of milk were mainly 
related to variation in aspects of farm performance and char¬ 
acteristics. Thus, this indicates that grass-based dairy 
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Table 4 Correlations (r) between various farm characteristics and carbon 
footprint (CF) of milk 


Farm characteristic 

CF of milk 

P value 

Farm size, ha 

0.07 

NS 

Cows, n 

-0.04 

NS 

Cows/ha 

-0.33 


Stocking rate, LU/ha 

-0.32 

*** 

Total farm milk production, t FPCM 

-0.17 

NS 

Annual FPCM yield/cow 

-0.48 


Annual FPCM yield/ha 

-0.55 

*** 

Fat, % 

-0.24 

** 

Protein, % 

-0.31 


Lactose, % 

-0.18 

NS 

Replacement rate, % 

0.03 

NS 

Grazing season, days 

-0.45 

*** 

Concentrate, kg DM/cow 

0.15 

NS 

Purchased forage, kg DM/cow 

0.20 

* 

Total feed, kg DM/cow 

-0.24 

** 

On-farm forage use, kg/ha 

-0.48 


On-farm N fertiliser, kg/ha 

0.09 

NS 

Farm-gate N surplus a , kg/ha 

0.28 

** 

Farm N efficiency 13 , % 

-0.66 

*** 

Electricity, kWh/cow 

0.10 

NS 

On-farm fuel, 1/ha 

0.06 

NS 


Fat and protein-corrected milk standardised to 4 % fat and 3.3 % tme 
protein per kilogram 

L U livestock unit (equivalent to the annual N excretion of an average Irish 
dairy cow), DM diy matter, NS not significant 

*P<0.05; **P< 0.01; ***P<0.001 

a N imports-N exports passing in or out through the farm-gate 
b Farm N exports/farm N imports 

producers can mitigate the CF of milk by adopting manage¬ 
ment practices that improve efficiency and performance. Fur¬ 
thermore, the study suggests that the goals of maintaining 
farm profitability and reducing CF of milk are not contradic¬ 
tory, given that previous studies of grass-based research farms 
report that improving farm productivity increases profit 
(Lovett et al. 2008; O’Brien et al. 2010). However, no single 
farm attribute accounted for the majority of variation between 
farms’ CF of milk. Thus, the study implies, similar to Beukes 
et al. (2010), that a suite of farm practices are required to 
increase efficiency and thus reduce carbon footprint of milk. 

4.1 Comparisons with national and international studies 

Relative to previous single-point estimates of the national 
average CF of Irish milk (Casey and Holden 2005; Lovett 
et al. 2008; Teagasc 2011), our mean result for commercial 
dairy farms was lower by 7-16 %. This was primarily because 
we considered a subset of dairy farms that were on average 



kg C0 9 -eq/kg FPCM 

Fig. 2 Histogram and expected normal distribution {bell-shaped line ) of 
carbon footprint of fat and protein-corrected milk (FPCM), weighted by 
farm FPCM production for 124 dairy farms in C0 2 equivalent (C0 2 -eq). 
The middle vertical solid line indicates the mean value and the dotted 
vertical lines to the left and right of the mean indicate the interval between 
the lower 2.5 % and upper 97.5 % of values. The mean carbon footprint of 
milk and standard deviation (SD) are also shown 

moderately more productive (e.g. higher milk yield/ha) than 
the national average farm. For example, compared to the 
performance of the average Irish dairy farm for 2010 and 
2011 (Teagasc 2011), the farms analysed produced more 
FPCM/cow (5,380 versus 5,058 kg) and per hectare (9,494 
versus 8,982 kg) and fed approximately 276 kg DM less 
concentrate/cow. In addition, the mean period cows spent at 
pasture was longer than the national average, which Lovett 
et al. (2008) reported and the present results indicate to miti¬ 
gate GHG emissions from grass-based dairy production. 
However, relative to a recent European estimate of the average 
CF of Irish milk by Leip et al. (2010), our results were 

Table 5 Regression coefficients (6), associated standard errors (SE), 
significance of regression coefficient and coefficient of determination 
(R ) estimated in a stepwise multiple regression model of carbon footprint 
of milk 


Parameter 

b 

SE 

P value 

R 2 

Intercept 

2.19 

0.08 

*** 

— 

N efficiency a , % 

-0.98 

0.10 

*** 

0.46 

Grazing season, days 

-2.16 x 1(T 3 

2.71 x 1(T 4 

** 

0.61 

FPCM yield/cow 

-7.76 x 10 -5 

9.87 xlO -6 

** 

0.70 

Replacement rate, % 

3.37 x 1(T 3 

9.43 xlO -4 

* 

0.75 


Fat and protein-corrected milk standardised to 4 % fat and 3.3 % true 
protein per kilogram 

*P<0.05; **P<0.01; ***P<0.001 

a Farm N exports/farm N imports 
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approximately 10 % higher, despite the greater efficiency of 
the dairy farms we assessed. The anomaly, though, was simply 
due to different LCA modelling assumptions. 

For instance, Leip et al. (2010) included carbon sequestra¬ 
tion by permanent grassland, which PAS 2050 excludes based 
on the IPCC (2006) guidelines. Recent studies report that 
permanent grasslands are an important long-term carbon sink 
(Soussana et al. 2007, 2010). Adopting the same carbon 
sequestration assumptions in the present study as Leip et al. 
(2010) reduced the mean CF of milk to 0.95 kg of C0 2 -eq/kg 
FPCM, which is 5 % lower than the estimate of Leip et al. 
(2010) of the average CF of Irish milk. Thus, as more carbon 
sequestration data becomes available, PAS 2050 may need to 
be revised to include this sink. The large influence that alter¬ 
native modelling choices (e.g. allocation methodologies) have 
on the CF of milk and meat is a well-documented LCA issue 
(Flysjo et al. 2011). Thus, until a harmonised LCA approach 
such as PAS 2050 is widely adopted internationally, direct 
comparisons between CF studies of milk are of limited value. 
Consequently, the FAO is leading a livestock environmental 
assessment and performance partnership (LEAP 2014) to 
develop global LCA guidelines for a wide array of environ¬ 
mental impacts from livestock. The draft LEAP guidelines 
generally adopt the same principles as PAS 2050, but are 
specific for animal supply chains. 

Although methodological variations often partially explain 
differences between LCA studies, cautious comparisons can 
provide an indication of the validity of LCA results and are 
useful in understanding the potential to mitigate CF of milk. 
Our results were similar to pan-European and global LCA 
studies, which showed that grazing systems of developed 
nations have the lowest CF of milk (Leip et al. 2010; 
Hagemann et al. 2012; Opio et al. 2013). However, such 
studies generally only assess the national or regional average 
situation and do not assess the variation in CF of milk within a 
region. Several reports suggest that this variation is 
significant. For instance, DairyCo (2012) showed that even 
though the average CF of UK milk was in the lower range of 
literature estimates, substantial differences were found 
amongst 415 grazing and non-grazing dairy systems’ CF of 
milk (0.83-2.81 kg of C0 2 -eq/l of fat corrected milk). A 
similar variability has also been reported for CF of Dutch milk 
(Thomassen et al. 2009) and USA milk (Thoma et al. 2013) 
where commercial non-grazing systems were predominately 
assessed and by an Australian evaluation of 140 grazing and 
non-grazing dairy systems (Dairy Australia 2012). In addition, 
the results of these studies show that farmers of developed 
nations can reduce CF of milk via management changes, e.g. 
DairyCo (2012) reported that reducing concentrate supple¬ 
mentation of cows decreased CF of milk of dairy farms. Thus, 
this indicates that global or European single-point estimates of 
national CF of milk may underestimate the scope to mitigate 
GHG emissions. 


Generally, our analysis agreed with large-scale LCA stud¬ 
ies (e.g. Thomassen et al. 2009; Dairy Australia 2012) regard¬ 
ing the range of farms’ CF of milk and the main sources of 
dairy farms’ GHG emissions (e.g. enteric CH 4 ). However, 
unlike LCA studies of similar scale, the farms we considered 
were not nationally representative and, as discussed margin¬ 
ally, more productive than the national average dairy farm. 
Consequently, the variability amongst farms’ CF of milk was 
lower than previous large-scale LCA studies (e.g. DairyCo 
2012). There was, however, a significant variation amongst 
key sources of dairy farms’ GHG emissions, particularly N 2 0 
emissions from manure and emissions associated with on- 
farm fertiliser use. The variability of these GHG sources was 
similar to previous reports by Dairy Australia (2012) and 
Thoma et al. (2013). Thus, these results suggest that there is 
potential to mitigate GHG emissions of grass-based dairy 
systems, which can be realised in part by adopting new 
technologies, but also via changes in farm management 
practices. 

4.2 Impact of farm efficiency and performance on CF of milk 

Similar to Casey and Holden (2005), Christie et al. (2011) and 
Yan et al. (2013b), numerous measures of farm efficiency and 
performance were correlated to CF of milk from grass-based 
farms. Congruous with Yan et al. (2013a), the farm attribute 
that had the strongest association with CF of milk was farm N 
efficiency (farm N imports/farm N exports). Key management 
practices that improve N efficiency of grass-based farms 
included adoption of white clover, greater utilisation of 
manure and better timing of manure and fertiliser application 
to grass growth. These practices reduced surplus N or the 
requirement for N fertiliser, which increased farm N 
efficiency. However, reducing N use on an area or animal 
basis rather than per unit of milk had little or no influence on 
CF of milk, which agrees with similar analysis by Olesen et al. 
(2006). Therefore, this may explain why several studies have 
generally not found or reported weak relationships between 
farm N use and CF of milk, given that most studies only assess 
farm N use per hectare (Dairy Australia 2012; DairyCo 2012). 

Grass-based farms that were not efficient from an N per¬ 
spective, however, did not necessarily have a high CF of milk. 
The main explanation for this was improving farm N efficien¬ 
cy primarily reduced emissions associated with N fertiliser, 
but had no effect on enteric CH 4 emissions, which were the 
main components of farms’ CF of milk. Therefore, farms that 
emitted lower enteric CH 4 emissions/unit of milk than more N 
efficient farms could achieve a lower CF of milk. Increasing 
milk yield/cow and reducing the annual herd replacement rate 
were the primary farm attributes that mitigated enteric CH 4 
emissions. Improving these attributes increased feed conver¬ 
sion efficiency (total feed DM/kg of FPCM), which in agree¬ 
ment with Rotz et al. (2010) was the main determinant of 


<£) Springer 



Int J Life Cycle Assess (2014) 19:1469-1481 


1479 


enteric CH 4 emissions/unit of milk. Increasing cow genetic 
merit via artificial insemination was the main management 
practice available to farmers to increase cow milk perfor¬ 
mance. Furthermore, the practice could also be used to im¬ 
prove cow fertility or health, which reduces the requirement 
for replacement heifers. However, this was achievable only 
using sires of a total genetic merit index, e.g. Irish economic 
breeding index (EBI), which are bred to increase cow perfor¬ 
mance, fertility and health. 

Improving milk yield/cow also reduced CF of milk, but 
only caused minor reductions in GHG emissions from N 
inputs, e.g. fertiliser. In addition, similar to Ramsbottom 
et al. (2012), higher milk yield/cow was associated with 
greater concentrate feeding, which increased GHG emissions. 
Consequently, milk yield/cow was not as influential in deter¬ 
mining CF of milk as farm N efficiency. Improving milk yield 
per hectare rather than per cow had less of an influence on feed 
conversion efficiency, but tended to cause a greater improve¬ 
ment in N efficiency. Furthermore, the measure increased 
homegrown forage utilisation per hectare, which in agreement 
with previous studies, reduced emissions from concentrate 
(Thomassen et al. 2009; O’Brien et al. 2010). Thus, 
improving milk yield per hectare had a slightly greater 
mitigating influence than increasing milk output/cow on CF 
of milk. Stocking rate was a key determinant of milk yield per 
hectare. However, stocking rate had little or no effect on feed 
conversion efficiency. Thus, this explained the lower influ¬ 
ence milk yield per hectare had on feed conversion efficiency 
relative to milk yield/cow and the weak negative association 
between stocking rate and CF of milk. 

Homegrown forage utilisation per hectare was primarily 
determined by stocking rate, but had a stronger association 
with CF of milk. This was because the farm attribute, unlike 
stocking rate, tended to be positively associated with the 
grazing season. Extending the grazing season shortened the 
winter housing period, which similar to Schils et al. (2005, 
2007) reduced GHG emissions from on-farm fuel use, manure 
storage and spreading. In addition, the management practice 
increased grazed grass utilisation per hectare, which decreased 
GHG emissions from concentrate and, in contrast to the 
results of Schils et al. (2007), reduced enteric CH 4 emissions. 
However, the main forage fed to cows indoors in the study of 
Schils et al. (2007) was maize silage, whereas in the present 
study, grass silage was mainly fed, which yields more enteric 
CH 4 than grazed grass (Robertson and Waghom 2002). There¬ 
fore, our results showed, unlike Schils et al. (2005,2007), that 
extending the grazing season reduced CF of milk, but the 
negative association was only moderate. This was because 
the farm practice also increased N 2 0 emissions from manure 
excreted on pasture by grazing cattle. 

No single farm attribute analysed explained the majority of 
variation between the CF of milk from grass-based farms. As a 
result, similar to previous studies, improving farm measures in 


isolation had a minor mitigating effect on CF of milk (Lovett 
et al. 2008; Dairy Australia 2012; Dairy Co 2012). Neverthe¬ 
less, CF of milk from grass-based farms could be well pre¬ 
dicted using just farm N efficiency, the length of the grazing 
season, milk yield/cow and the annual herd replacement rate. 
Thus, it may not be necessary to conduct detailed on-farm 
surveys to estimate CF of milk. However, to validate the 
prediction capabilities of these farm attributes, a new indepen¬ 
dent audit would need to be conducted. 

As discussed, the farm attributes that predicted CF of milk 
were influenced by farm management practices, e.g. stocking 
rate, which varied to a similar or greater extent as GHG 
emissions between farms. Therefore, in order to reduce carbon 
footprint of milk, grass-based dairy farmers need to implement 
a suite of management practices to simultaneously improve 
farm N efficiency, the length of the grazing season, milk yield/ 
cow and the annual herd replacement rate. Farm management 
practices that positively influenced these key determinants of 
CF of milk and did not negatively interact were improving 
cow total genetic merit via artificial insemination, extending 
the length of the grazing season, reducing concentrate feeding, 
increasing stocking rate and reducing fertiliser N use per 
hectare by increasing the proportion of manure applied in 
spring. Furthermore, implementing these practices improves 
farm productivity. This increases the economic viability of 
grass-based farms (Lovett et al. 2008; Beukes et al. 2010; 
O’Brien et al. 2014) which is a key consideration when 
changing management practices. Therefore, this implies that 
grass-based dairy farmers can potentially significantly miti¬ 
gate CF of milk and maintain farm profitability. 


5 Conclusions 

Independent certification of an LCA model according to PAS 
2050 for a large group of grass-based dairy farms provides a 
verifiable approach to quantify CF of milk at a farm or 
national level. The application of the certified model in this 
study showed that the mean CF of milk from a subset of grass- 
based Irish farms was in the lower range of literature 
estimates. However, comparisons with previous LCA studies 
that were not certified to the PAS 2050 standard were affected 
by inconsistent modelling assumptions and choices. 
Therefore, this highlights the need for experts to use an 
internationally standardised LCA approach to estimate CF of 
milk, for instance PAS 2050 or IDF (2010) dairy LCA guide¬ 
lines, which can be externally certified by an accredited third 
party, e.g. Carbon Trust. 

As expected, differences between grass-based farms’ CF of 
milk were primarily explained by variation in measures of 
farm performance, particularly farm N efficiency. Therefore, 
this suggests grass-based farmers can reduce CF of milk and 
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maintain or increase farm profit by adopting management 
practices that increase efficiency, e.g. improving cow total 
genetic merit. However, no individual farm attribute analysed 
explained the majority of variation between the CF of milk 
from grass-based farms. Therefore, the study indicates that 
several farm practices are required to reduce CF of milk from 
grass-based farms. 
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